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ABSTRACT 

In several recent observational studies on Type Ib/c supernovae (SNe Ib/c), the inferred 
ejecta masses have a peak value of 2.0 - 4.0 Mq, in favor of the binary scenario for their pro¬ 
genitors rather than the Wolf-Rayet star scenario. To investigate the observational properties 
of relatively low-mass helium stars in binary systems as SN Ib/c progenitors, we constructed 
atmospheric models with the non-LTE radiative transfer code CMFGEN, using binary star evo¬ 
lution models. We find that these helium stars can be characterized by relatively narrow helium 
emission lines if the mass-loss rate during the final evolutionary phase is significantly enhanced 
as implied by many SN Ib/c observations. The optical brightness of helium star progenitors can 
be meaningfully enhanced with a strong wind for M > 4.4 Mq, but hardly affected or slightly 
weakened for relatively low-mass of ^ 3.0 Mq, compared to the simple estimate using blackbody 
approximation. We further confirm the previous suggestion that the optical brightness would be 
generally higher for a less massive SN Ib/c progenitor. In good agreement with previous studies, 
our results indicate that the optical magnitudes and colors of the recently detected progenitor of 
the SN Ib iPTFlSbvn can be well explained by a binary progenitor with a final helium star mass 
of about 3.0 - 4.4 Mq. 

Subject headings: binaries: general — stars: atmospheres — stars: evolution — supernovae: general — 
supernovae: individual (iPTFlSbvn) 


1. Introduction 

It is well known that the majority of massive 
stars are born in multiple systems, and binary 
systems play an important role for core-collapse 
supernovae (SNe). Recent observations indicate 
that a large fraction of massive binary stars are 
in relatively short-period orbits, implying that 
they would undergo binary i nteractions durin g 
the course of their evolution ( Sana et al ][20 i 3). 
Furthermore, the mass ratios of the stellar com¬ 
ponents of massive binary systems are found to 
be close to one for many cases, providing a favor¬ 
able con dition for stable mass tra nsfer to avoid 
mergers ( Kobnlnickv fc Frverll200^ . These flnd- 
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ings strengthen the long-standing argument that 
binary interaction should be considered one of the 
primary factors for massive star evolution, and 
that binary stars may not only be related to cer¬ 
tain exotic SNe, but also to commonly observed 


ones ii 

rcluding Tvpe Ib/c 

supernovae (SNe Ib/c) 

(e.g. 


Posdiadlowski et al. 

|l992|: IWooslev et al. 

1995 


Wellstein & Langeii 

1999: Eldridge et al. 

2008 

:lYoon et al. 201f)l Smith et al. 201lll . 


SNe Ib/c, whi ch constitute about 25% of all 
core- collapse SNe ( Smith et al.l20Il ; Eldridge et al 


20131) . are therefore important to constrain the 
evolution of massive stars. The hydrogen en¬ 
velopes of the progenitors of these SNe must have 
been stripped off before they exploded. Both 
mass loss from single stars and binary interac¬ 
tions can fulfill this condition, but the result¬ 
ing properties of SN Ib/c progenitors would be 
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very different from e ach other ( Yoon et al ] l2012l: 
Eldridge et al ] |2013l) . Single Wolf-Rayet (WR) 


stars in the nearby Uniyerse haye bolometric lu¬ 
minosities of \ogL/LQ ^5.1, implying that their 
masses are higher than about 10 Mq. Their final 
masses at the pre-SN stage are predicted to be 
higher than about 7 Mq. By contrast, lower final 
masses are strongly preferred in binary progeni¬ 
tors. Recent analyses on the light curyes of SNe 
Ib/c indicate that SN Ib/ c ejecta masses are typi¬ 


cally around 1.0-5.0 Mq (|Drout_et_^jJl20n ; 
120131 : Ibyman et al. 20141 iTaddia et al.1l2015 ) 


Canol 


sup¬ 
porting the binary scenario oyer the WR star sce¬ 
nario. 

Howeyer, the debate on which type of pro¬ 
genitors between single WR stars and relatiyely 
low-mass helium stars in binary systems is the 
dominant o ne for SNe Ib/c is still on-going (see 
Yoon 20151 for a recent reyiew). The best way 
to resolye the issue would be therefore to di¬ 
rectly identify S|£^Ib/c progenitors in pre-SN im¬ 
ages (Smartt 2009ll . The previous searches haye 


not been successful (Crockett et al.ll2007t ISmartt 


20091 Eldridge et al. 2013l l. except for the tenta¬ 


tive identific ation of the pro genitor of the SN lb 
iPTFlSbvn (ICao et al.ll2^ . 

Yoon et ^ ( 2012ll made rough estimates on the 


optical magnitudes of SN Ib/c progenitors at the 
pre-SN stage using their evolutionary models of 
massive stars. They concluded that WR progen¬ 
itors would be generally very faint in the optical 
(My ss —2 ^ —3 mag) compared to most WR 
stars observed in the nearby Universe, while rela¬ 
tively low-mass helium star progenitors in binary 
systems can be much brighter in the optical, as 


they b ecome helium giant stars. Eldridge et al.l 


(|2015l) also made a similar conclusion on binary 
progenitors. But this conclusion is based on the 
simple assumption of blackbody radiation from 
the adopted stellar evolution models. In reality, 
both absorption and emission lines from helium 
stars may have a significant impact on the optical 
brightness, depending on the surface temperature 
and the mass-loss rate by winds. 

( 2013b[) and Groh et ^ ( 2013a 1 


Groh et al 


presented stellar atmospheric models of single WR 
type progenitors at the pre-SN stage using the stel¬ 
lar evolution models of the Geneva group. They 
concluded that WO type progenitors, of which 
the initial masses are higher than about 35 Mq, 


would be relati vely faint in the o ptical, in good 


agreement with lYoon et al.l (j2012n . and that WN 


type progenitors, which are expected for a lim- 
ited initial mass ran ge (M = 31 — 35 Mq; but see 


Eldridge et al.ll2015ll . can have optical magnitudes 


of about — 5.5 ^ —6.5 mag. However, there exists 
no such a study on binary progenitors. 


The purpose of this paper is therefore to present 
stellar atmospheric models of binary progenitors 
of SNe Ib/c, for the first time. This would al¬ 
low to better test the binary progenitor scenario 
with iPTFlSbvn as well as future observations. In 
Sect.O we explain our numerical methods and the 
adopted physical assumptions. The result of our 
calculations is summarized in Sect. |31 We com¬ 
pare our model spectra with those of HD 45166 
and V Sgr which are relatively low-mass helium 
stars observed in our Galaxy, as well as WN stars 
in Sect. IH The predicted optical magnitudes and 
the implications for the progenitor of iPTF13bvn 
are discussed in Sects. [5] and IH We conclude this 
study in Sect. [T] 


2. Physical assumptions 


The helium star progenitor models for the 
present study were chosen f rom t he binary evolu¬ 
tion models bv lYoon et al.l ( 20in[l . These models 
were calculated at solar metallicity including the 
effects of rotation, and followed up to the end of 
core neon burning, which is about one year be¬ 
fore core collapse. The surface properties of these 
models remain almost unchanged during the fi¬ 
nal evolutionary stages after carbon exhaustion, 
which last for about 100 years. The pre-SN images 
used for the searches of SN Ib/c progenitors were 
usually ta ken about 10 years be f ore the SN explo - 
sion (e.g., Grockett et al ] 120071: iGao et al.l lirTil 
and therefore any model after core carbon exhaus¬ 
tion are suitable for comparison with observations. 
The models chosen for atmospheric calculations 
are the last computed models of the sequences 
5, 22 and 27, which have the final masses of 3.0, 
4.4, and 5.1 Mq, respectively. The initial masses 
of these models are 14, 18, and 25 Mq, respec¬ 
tively. The optical luminosity of a helium star at 
the pre-SN stage is a sensitive f unction of its mass 


and radius ( Yoon et'al] 2012 ). Yoon et al. ( 20I0h 


showed that there exists a fairly good mass-radius 
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relation for binary SN Ib/c progenitors at the 
pre-SN stage, and therefore these selected models 
can roughly represent the surface properties of SN 
Ib/c progenitors for the given final masses. 

Given the neutron star remnant mass of about 
1.3 Mq, this range of helium star masses is con¬ 
sistent with the peak values of the SN Ib/c ejecta 
mass distribution (i. e., Mf.\ fti 2 — 4 inferred 
from SN light curves (IProut et aLll201ll : ICao et al 


20131 : iLvman et al.l l2014ll . It also fits well with 
the infe rred ejecta mass of iPTF13bvn {Mf.\ s a 
2.0 Mq-, Fremling et al.l[2014 : Bersten et al.ll20M) . 
and thus we can directly compare our result to the 
observed properties of the tentative iPTF13bvn 
progenitor. 


Yoon et alJ (120101) adopted the WR mass-loss 
rate given bv iHamann et ah (1995) using a reduc¬ 


tion factor of 5 - 10 to consider the effect of WR 
wind clumping. However, the resulting wind mass- 
loss rates of these models are not high enough to 
have an optically thick WR type wind, in partic¬ 
ular for 3.0 and 4.4 Mq helium stars, given that 
their masses are relatively low and that their en¬ 
velopes have relatively large radii. But the mass- 
loss rates from helium stars of the considered mass 
range are not well constrained by observations. 
These helium stars undergo very rapid increase 
in the surface luminosity as they ap proach to the 
pre-SN stage ( Yoon et'^l2010ll201^ . and we can¬ 
not exclude the possibility that mass loss becomes 
dramatically strong during this final stage (i.e., 
M > 10“® Mq yr“^), which is in fact implied 


by inany SN Ib/c observati o ns le.g.. iFoley et al 
2007: ^storello et al.1 12 OO 8 : Wellons et al.l 20ia 


Gorbikov et al ■11201 4^ . Therefore, we also consider 
an arbitrarily enhanced mass-loss rate for each he¬ 
lium star model as summarized in Table [TJ Here, 
the label wl den otes the models w ith the mass- 
loss rates used bv lYoon et al. ( 2010ll . and w2 the 
models with the enhanced mass-loss rates. 

The spectra of our SN Ib/c progenitor mod¬ 
els were computed using the non-local thermody¬ 
namical equilibrium atmospheric radiative trans- 
fer code CMFGEN ( Hillier fc Milleil|l998l : lllillieij 
2003). CMFGEN determines the temperature dis¬ 
tribution of the expanding atmosphere by solving 
the statistical and radiative equilibrium equations, 
and computes line and continuum f ormation wit h 
the spherical symmetric geometry ( Hillieil 1990(l . 
CMFGEN uses the super-level approach to fully 


treat line blanketing. In this approach, levels with 
similar excitation energies are grouped into a sin¬ 
gle level, under the assumption that the departure 
coefficients in a group are identical, and only the 
population of the super level is solved to to specify 
the populations of the levels withi n a super level 
( Hillier fc Miller 19981: Hillieil 20031) . 

CMFGEN requires a complete (previously con¬ 
verged) model including atmospheric structure, 
atomic models and their departure coefficients as 
the initial trial solution. For our calculations, we 
first adopted one of the pre-calculated O star mod¬ 
els provided by Hilliei0 with the effective tempera¬ 
ture of 27,500 K and Log g of 3.25, which are sim¬ 
ilar to those of our 4.4 Mq helium star (Table [TJ. 
Table [2] presents the atomic species included in our 
calculations with the numbers of super and full 
levels. The atomic data and information on the 
levels are accompanied with the CMFGEN code. 
Once we obtained the converged model for 4.4 Mq 
helium stars, we used it as the initial trial solution 
for the other models. The atomic models from 
the above O star do not include neutral and low- 
excitation level metal lines, which may be relevant 
for our 3.0 Mq model, but this does not signifi¬ 
cantly affect the main conclusions of our study as 
explained in Section [4] below. 

While a hydrostatic solution for the subsonic 
part is self-consistently solved by CMFGEN, it 
does not solve the momentum equation of the 
wind, and thus the wind velocity profile needs to 
be specified. For the wind part, we assume the 
standard /3-type velocity law. We used /3 = 1 and 
1.5 and Uoo/'Uesc = 2.0 and 1.5 for the optically thin 
wind models (wl) and t he WR-type wind mod - 
els (w2), respectively ('cf. IVink fc de Koterll2005 b 
where Voo and Vesc respectively denote the termi¬ 
nal wind velocity and the escape velocity. The ve¬ 
locity structure is modified at depth to smoothly 
match the structure at the surface of the hy¬ 
drostatic core (for details on the CMFGEN cal 


culations, see Hillier 19901 : Hillier fc Miller Il998t 


Hillieil 1200311 ■ 


Yoon et al. ( 2010ll considered massive compan¬ 
ion stars of OB type in their calculation. The 
companion star masses of the chosen progenitor 
models (M = 3.0,4.4, and 5.1 Mq) at the end 
of the calculation are 18, 23, and 35 Mq, re- 


^ http: //kookaburra.phyast .pitt .edu/hillier /web / CMFGEN. htm 
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specti vely. Note that th ese are not unique solu¬ 
tions. Yoon et al. ( 20inll followed the spin-up of 
mass accreting star as a result of mass and an¬ 
gular momentum accretion and the resulting en¬ 
hancement of mass loss, and thus the mass transfer 
process is highly non-conservative in their calcu¬ 
lations. However, the companion star masses can 
significantly vary according to the adopted mass 
accretion efficiency (i.e., the ratio of the trans¬ 
ferred mass to the accreted mass). For exam¬ 
ple, if the mass transfer were conservative, a bi¬ 
nary system where both stars have an initial mass 
of about 20 M© could produce a 3 M© helium 
star that has a structure similar to that of the 
3.0 Mq models of the present study, via Case A 
and AB mass transfers. The corresponding com¬ 
panion s tar mass would be about 35 M© in this 
case Icf. IWellstein fc LangeilllO^ . 

Given that many different combinations of he¬ 
lium star and companion star masses are possible 
in principle, we included 0-type star models of 
three different masses (M = 20, 30, and 35 Mq) 
in the CMFGEN computation to investigate how a 
luminous 0-type star companion would influence 
the optical brightness of SN Ib/c progenitor sys¬ 
tems. These masses were chosen mainly because 
they can give the best fit with the observational 
properties of the iPTF13bvn progenitor candidate 
as discussed below. The 20 and 30 M© models 
are non-rotating models without overshooting in 
their early stages of the main sequence that are 
close to the zero-age main sequence. The 35 M© 
model is the last com puted companion s tar model 
of the sequence 27 of Yoon et al. ( 2010l l (i.e, the 
sequence for the 5.1 M© helium star model in the 
present study). 

The model parameters are summarized in Ta¬ 
ble □ The table lists two kinds of temperature 
and radius. The values of T* an d Rj, are predicted 


from the stellar evolution code ( Yoon et al. 20inl l 


without correcting the optical depth effects from 
the wind, whereas the effective temperature (Teff) 
and photospheric radius (i?phot) are the outputs 
from the GMFGEN calculations defined at the 
layer where the Rosseland optical depth (rross) is 
2/3. In the calculations, hydrodynamic clumping 
in the wind is considered by adopting a volume fill¬ 
ing factor (/). For all models, we assume / = 0.1. 
The surface abundances were taken from the data 
of the selected evolutionary models. 



Fig. 1.— SEDs of the helium star models and 
the blackbody fluxes (purple dashed line) of the 
given temperature T* and luminosity of each mass 
model. All fluxes are scaled to a distance of 10 pc. 
Transmission curves of the HST/ACS Wide Field 
Ghannel filters, F435W, F555W, and F 8 I 4 W are 
also presented at the bottom panel. 


3. Results of the atmospheric models 

Spectral energy distributions (SEDs) of the he¬ 
lium stars from the GMFGEN calculations are 
shown in Figure [1] All fluxes are scaled to a 
distance of 10 pc. The SEDs of the wl and w2 
cases for a given helium star mass are almost sim¬ 
ilar to each other, but the w2 models have ex¬ 
cess in long wavelengths because of the extended 
photosphere radius with an optically thick wind. 
For comparison, we also present the blackbody 
fluxes of the given temperature (T*) and luminos¬ 
ity for each model. In optical wavelength ranges, 
the differences between the blackbody fluxes and 
the helium star models are not significant. How¬ 
ever, the w2 models, particularly for M = 4.4 and 
5.1 Mq, show strong excess in continuum as well 
as a wealth of emission lines in the optical. The 
large differences in the extreme ultraviolet wave- 
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Fig. 2.— Normalized optical spectra of the helium 
star and 0-type star models for the wavelength 
range of 3750 A < A < 4900 A. 
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Fig. 3.— Normalized optical spectra of the helium 
star and 0-type star models for the wavelength 
range of 5750 A < A < 7200 A. 


lengths are owing to the line blanketing effect. 

Figures [2] and [3] present the normalized opti¬ 
cal spectra of all models. While the spectra of 
0-type stars only show absorption lines, those of 
helium stars (particularly the w2 models) show 
emission lines in general. The 3.0 Mq models 
and the 4.4 Mq model with wl type wind (M3wl, 
M3w2, and M4wl) dominantly show absorption 
lines; however, compared to 0-type stars. He i 
lines of the helium stars are broader with devel¬ 
oped wings (e.g., He i A4024 or He i A4472), and 
they sometimes appear as (weak) P Cygni pro¬ 
files (e.g., He I A6679 or He i A7067). In addition, 
the Ha line of the M4wl model appears in emis- 
sioiil. These emission lines make the helium star 
to be distinguishable from O stars with a similar 
temperature and mass-loss rate. The H lines in 
the M3wl model are absorption lines and weaker 
than those of 0-type stars because of the lower 
abundance of hydrogen (see Table [T]). The spec¬ 
tra of the 5.1 Mq models of which hydrogen is 

^ As discussed by lYoon et al.l 1I2OIOI I in detail, the sur¬ 
face mass fraction of hydrogen for the 3.0 and 4.4 Mq 
models is about 0.01 (Table [T}. Hydrogen lines are 
actu ally detected in many SNe classified as Type Ib 

(e.g., [peng et al.|200Cl|;|Branc h et al.l2002l : lElmhamdi et al.l 

120061 : ISpenceTfc^ar^ 120101 1! and this theoretical pre¬ 

diction of the tiny amounts of hydrogen in some SN 
Ib/c progenitors is consistent with observations (see also 
iDessart et al.|[201lf) . 


almost depleted show broad helium emission lines 
and high-ionization lines such as N iii/iv. 


4. Comparison with Observational Coun¬ 
terparts 

One of the observational counterparts of binary 
SN Ib/c progenitors is the quasi-WR (qWR) star 
HD 45166, which is composed of a helium rich 
4.2 Mq star with i? ~ 1.0 Rq and a 4.8 Mq 
main sequence star (B7V) in a 1.596 day orbit 


(Willis fc Sticklan3l983: Steiner &: Oliveir J2005 


Groh et al.l 120081) . The optical spectrum of the 
helium star that is likely on the helium main se¬ 
quence shows a number of emission lines such as 
He i/ii, N iii/iv/v, and C m/iv, and it is well ex¬ 
plained by a model with Teff = 50, 000 ± 2, 000 K 
(T* = 70,0p0± 20,000 K), logfL/L ^) = 3.75 ± 
0.08, and M = 2.2 x 10“^ Mq yr“^ ( Groh et al. 
I 2 OO 8 I ). The spectrum of the qWR star in HD 45166 
with weak He i lines and higher-ionization lines 
(e.g., N IV or C iv) is different from our helium star 
model spectra of similar mass (M4wl and M4w2), 
which dominantly show He i lines. This difference 
is likely due to the much lower effective tempera¬ 
ture (Teff = 20,000- 30,000 K) o f our 4 . 4 Mq mod¬ 
els. As discussed in Yoon et al. ( 2010l . 2012 ). low- 
mass (M = 3-5 Mq) helium stars are hot and vi¬ 
sually very faint on the helium main sequence but 
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become cooler and more luminous during the final 
evolutionary stages because of the rapidly expand¬ 
ing envelopes. The qWR star in HD 45166 is in 
fact much fainter in the o ptical {My = —0.21 mag; 
IWillis fc Stickland 1983ll than our 4.4 Mq models 
{My fv —5 mag, see Section 3). 

Another observational counterparts are evolved 
helium giant stars in binary systems. Only four 
stars of such a system are currently known in- 
cluding V Sgr, KS Per, LSS 1922, and LSS 4300 
(IDudlev fc JeffervI 1199311 . and the v Sgr system 
among them has been best studied. The primary 
star of the v Sgr system is a hydrogen-deficient 
star with M ^ 3.0 Mc^, T^ff ^ 12,000 K, and 
logL/ ~ 4.6 ( Said 1995 : Kipper k. Klochkova 


20121 ) . which are comparable to those of our 


3.0 Mq helium star models. The main character¬ 
istics of the V Sgr spectrum includes He i absorp¬ 
tion lines, a large number of absorption/emission 
lines of neutral and ionized metals from low ex¬ 
citation levels, some P Cygni profiles, and the 
forbid den lines of low-ionization m etal lines (Fig¬ 
ure m iKipper fc Klochkoval l2012l ). He i absorp¬ 
tion lines and metallic lines such as C ii, N ii, 
or Fe III present similar characteristics for the 
V Sgr spectrum and the M3wl model spectrum. 
iBonneau et al.l ( 201 ill argue that the hydrogen P 
Cygni profiles of the v Sgr originate either from 
the circumbinary disk of this system or from the 
disk around the unseen secondary star. This im¬ 
plies that the wind mass-loss rate from the v Sgr is 
comparable to the value adopted for M3wl model 
that do not show emission lines, rather than that 
of M3w2 model for which the spectrum is charac¬ 
terized by strong emission lines for both hydrogen 
and helium (Figures [5] and [3]). A number of other 
absorption lines only seen in the v Sgr spectrum 
are neutral and metal lines of low excitation levels 
(e.g., N I, Ne i, or Fe ii), which are not included 
in our models. 

Finally, the spectrum of the M5w2 model may 
be compared with the spectra of WN stars with 
similar temperatures. In the M5w2 model spec¬ 
trum (Figures [3] and [3]), helium and nitrogen lines 
are dominantly observed, and most helium lines 
are broad and appear as P Cygni profiles. These 
spectral characteristics can also be found i n the 
spectra of WN stars ( Hamann et al. 200^. In 
particular, WN8 stars with T* ~ 40,000-50,000 K 
show very similar spectra with strong He i lines 



Fig. 4.— Comparison of the v Sgr spectrum 
(black) that was obtained by a 60-centimetre tele¬ 
scope at the Observatoire du Pic du Midi on 
July 24, 2008, which is available from the BeSS 
database (http://basebe.obspm.fr/basebe), with 
the 3.0 Mq helium star model spectra (dark green 
and light green for M3wl and M3w2, respectively). 
Vertical lines at the top of each panel mark spec¬ 
tral lines seen in the both v Sgr and helium star 
models. Black, red, brown, blue, and purple lines 
correspond to H i. He i, C ii, N ii, and Fe ill, 
respectively. In the lower left panel, the M3w2 
spectrum is scaled to 50% for display. 


and negligibly weak C iv line (at 5803A) as seen 
in the M5w2 model spectrum. As moving to¬ 
ward earlier type WN stars (i.e., higher surface 
temperature). He i lines become weaker and the 
C IV A5803 line becomes stronger. Compared 
to WN8 s tar spectra (e-K-. WR 012, WR040, or 
WR170 in iHamann et al.ll2006f) . however, in the 
M5w2 model spectrum He i lines tend to be 
stronger, and the N ill A4635/4I line is stronger 
than the He ii A4687 line. In addition, in spite of 
similar temperature, the M5w2 model is fainter 
in the optical {My > —5 mag, see Section 3) 


than WN8 stars {My ~ —7 mag; iHamann et al 


20061 ). which is a natural consequence of the lower 


bolometric luminosity for the given surface tem¬ 
perature. Compared to the case of M5w2, in the 
spectrum of M5wl the lines from higher excita¬ 
tion species such as N iv or C iv are stronger and 
He I lines almost disappear likely due to higher ef¬ 
fective temperature compared to the M5w2 model 
(Table [T|). 
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Fig. 5.— Predicted absolute magnitudes of the models in the HST/ACS F435W (left), F555W (middle), 
and F 8 I 4 W (right) filters. Filled cyan diamonds without borders are the absolute magnitudes of the helium 
star models for given HST/ACS filters, and open diamonds of brown, yellow, and purple are the magnitudes 
of 20, 30, and 35 M© 0-type stars, respectively. The filled cyan diamonds with borders are the predictions of 
the helium stars in binary systems with an 0-type star companion of corresponding colors. The extinction- 


corre cted magnitude ranges of the iP TF13bvn progenitor candidate by the low {E{B—V) = 0.0715. ICao et al 


2013 ) and high {E{B — V) = 0.2147. IPersten et al. 2014 ) extinction values are presented with red and dark 


green dashe d boxes, respective ly. Gray dashed lines present the prediction of 32 Mq WR star progenitor 
model from lOroh et al.l ( 2013al) . 


5. Absolute Magnitudes of Helium Star 
Binary Progenitors 

In Table [3] and Figure [H we present the 
predicted magnitudes of the helium star mod¬ 
els for the HST/ACS filters of F435W, F555W 
and F 8 I 4 W, which roughly correspond to John¬ 
son B, Johnson V, and B road I filters, respec¬ 
tively ( Sirianni et al.l 12005 ). We computed the 
absolute magnitudes of the models in the VEGA- 
MAG system using the synthetic model spectra 
and the transmission curves of the HST/ACS fil¬ 
ters presented in Figure [TJ For a filter P with 
a transmission curve of P(A), the absolute mag¬ 
nitudes in the VEGAMAG system are given by 
VEGAMAG(P) = 

-2.51ogio[/P(A)EA(A)AdA/ / P(A)FA,vega(A)A dA], 
where Ex and pA.Vega are the model flux from 
the CMFGEN calculations and the flux of Veg a 
scaled to a distance of 10 pc ( Sirianni et al.l 2005 ). 
The reference spectrum of Vega is adopted from 
the SYNPHOT package distributed as a part of 
STSDA^. We computed the bolometric mag- 


^ http: / / www.stsci.edu / institute / software.hardware / stsdas / synphot 


nitudes (Mboi) as welL assu ming that the solar 
Mboi is 4.74 mag ( Cox 2000l) . and present bolo¬ 
metric corrections in a given filter P (BCp = 
Mboi - VEGAMAG(P)) in Table H 


For comparison, the magnitudes and bolometric 
corrections of the blackbody models are also pre¬ 
sented in Table [31 The deviation from the black- 
body model prediction depends on the strength 
of various lines, as well as the location of the ac¬ 
tual photosphere for a given wind mass-loss rate. 
For example, if the photosphere were significantly 
lifted up with an optically thick wind that pro¬ 
duces strong emission lines, the resulting luminos¬ 
ity in the optical would be much higher than in 
the corresponding case of blackbody. Otherwise, 
numerous absorption lines tend to decrease the vi¬ 
sual brightness compared to what the blackbody 
models predict. 


In general, for a given mass-loss rate, the opti¬ 
cal thickness of a stellar wind becomes large r for a 
smaller radius of the star (e.g., Langeij[i989l) . This 
explains why the difference between the stellar 
radius (i?*) and the photospheric radius (i?phot) 
becomes larger with a higher mass for a given 
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wind mass-loss rate as shown in Table |3l This 
also tends to make emission lines stronger. There¬ 
fore, the visual brightness of helium stars can be 
strongly influenced by winds, and the assumption 
of the blackbody may significantly underestimate 
the brightness of a helium star, as shown with 
M4w2 and M5w2 models. As also expected, for 
a given helium star mass, such difference becomes 
larger with a higher mass-loss rate: the w2 model 
tends to be brighter by up to mag than the wl 
model. 


On the other hand, M3wl, M3w2 and M4wl 
models have numerous absorption lines in their 
spectra, and their photosphere radii are not much 
different from the stellar radii (see Section [3] and 
Figures [2] and Is]). Their visual brightness is com¬ 
parable or a little fainter (~0.3 mag) than the 
corresponding blackbody models, because of the 
presence of absorption lines and the lack of emis¬ 
sion lines. The 3.0 Mq model could be some¬ 
what fainter than predicted in Table 0 if we in¬ 
cluded low-excitation species, which can produce 
more absorption lines in the spectra (see the above 
discussion on v Sgr in Section a. However, for 
more massive models, our result would not be 
affected with the inclusion of those species be¬ 
cause the surface temperature is too high for the 
low-excit ation lines to make an y significant im¬ 
pact (e.g. Gray fc Corballvll2009ll . 

Note that the optical brightness becomes sys¬ 
tematically lower for a higher helium star mass, 
and that the considered helium stars in the present 
study are predicted to be much brighter than 
massive single WR stars of WO type (Mzams ^ 
30 — 35Mq) at the pre-SN stage that would have 
Mfsvr — —3 mag |Yoon et 'aL 2012t iGroh et ^ 
2013bll . The deepest absolute BVR magnitude 


limits of SN Ib/c progenitors in pre-explosion 
i mages to date are between —4 and —5 mag 
( Grockett et al ] 20071 : lEldridge et al.l l iolill . and 
they are comparable to the faintest one in our 
helium star mod els M5wl. Our work confirms 
the conclusion bv lYoon et al. ( 2012 ) that the non¬ 
detection of SN Ib/c progenitors does not neces¬ 
sarily imply binary progenitors rather than mas¬ 
sive WR progenitors. To the contrary, binary pro¬ 
genitors are easier to be found in the optical than 
massive WR progenitors of WO type that would 
be the most com mon type of SN Ib/c progenitors 


from single stars (|Yoon et al.ll20I2t lEldridge et al 


I 2 OI 5 I) . The probability of detecting binary SN 
Ib/c progenitors can further increase with a lumi¬ 
nous 0-type star companion as shown in Figure[5l 
while a significant fraction of binary progenitors 
would have a less luminous dwarf star or a com¬ 
pact object as a companion ( Yoon 20I5h . 


6 . Comparison -with the Progenitor Can¬ 
didate of the Supernova iPTFlSbvn 


Cap et al.l (|2013l) have recently reported a ten¬ 


tative identification of the progenitor of the SN 
Ib iPTFlSbvn exploded in NGC 5806 (22.5 Mpc) 
from the pre-explosion HST/ACS images. The 
observed magnitudes of the object in HST/ACS 
F4-35W, F555W, and F81/W images range from 
26 to 27 ma g dependirig on the adopted photome - 
try method ( Cao et 311120131 : lEldridge et al.ll20I5l) . 
The inferred absolute magnitudes range from -5.3 
to -7.3 mag, given the unc ertainties of the e xtinc- 
tion towards NGC 5806. ICao et all (2013) mea¬ 
sured the Milky-Way (foreground) and the host 
galaxy reddening from the observed Na i D lines 
as E{B — R)mw = 0.0278 and E(B — R)host = 


0.0437, whereas iBersten et al. (2 0141) adopted 


E(B — R)mw = 0.0447 (jSchlaflv fc Finkbednei 


2 OIII) and derived if (R — R)host = 0.17±0.03 using 


an intrinsic-color law from a sample of observed 
SNe. In this paper, we adopt the photometry re¬ 
sults of lEldridge et al. ( 2015 ) and consider both 
values of reddening to compare the progenitor can¬ 
didate of iPTF13bvn with our progenitor models. 
In FigurelH we present the magnitude range of the 
iPTF13bvn progenitor candidat e dereddened by 
the low {E{B — V) = 0.0715 from ICao et al.l 2013) 
and h igh {E{B — V)= 0.2147 from iBersten et ^ 
2 OI 4 I) extinction values. 


As shown in the figure, if we ignore the contri¬ 
bution from a companion star, only the 3.0 Mq 
progenitor can marginally satisfy the observed 
magnitude range of the iPTF13bvn progenitor 
with the low-extinction. For the 4.4 Mq and 
5.1 Mq progenitor, the companion star should be 
more massive than 20 Mq and 30 M 0 , respec¬ 
tively. With a 35 Mq companion, all progenitor 
models give optical magnitudes compatible with 
those in the high-extinction case, but they are a 
bit too bright to ht with the low-extinction case. 
Therefore, for the considered mass range of he¬ 
lium stars, 35 Mq roughly gives the upper limit of 











































































the companion star mass. From the evolutionary 
point of view, a more massive helium star pro¬ 
genitor would systematically have a more massive 
companion if the binary system underwent stable 
mass transfers. For example, with stable Case B 
mass transfer systems 0, the upper limit of the 
companion star mass would be roughly about 27, 
35 and 38 Mq, for 3.0, 4.4 a nd 5.1 Mq helium star 
progenitors lsee lYoonll2015L for more details on the 
binary progenitor evolution). With Case AB mass 
transfer system^, this limit would be higher b y 
several solar masses Icf. lWellstein fc LangeiJl99^ . 
The lower limit of the companion star mass is zero: 
a helium star progenitor of the considered mass 
range would not have any companion if it were pro¬ 
duced via Case A mass transfer that leads to re¬ 


versal of the SN order as discussed bv iPolsI (|l994ll . 

Considering both the theoretical constraint and 
the rather large magnitude uncertainty (~1 mag) 
resulting from the two extinction values, we con¬ 
clude that a helium star progenitor of 3.0/4.4 Mq 
with a 20/30 Mq 0-type star companion can best 
explain the observed magnitudes of the progenitor 
candidate in all the three HST/ACS filters. Their 
(F435VF-F555W) and {F555W-F8UW) colors 
presented in Figure [5] are also within the observed 
color ranges of the progenitor candidate, although 
all of the helium star binary progenitor models are 
distributed in a narrow color range. 

Groh et al. ( 2013all suggested a single WR star 


of WN type with initial masses of 31-35 Mq as a 
progenitor of iPTF13bvn based on the predicted 
optical magnitudes (My ^ —5.5 mag, gray dashed 
lines in Figure [5]) of the non-rotating models from 
the Geneva stellar evolution code. We note that 
Groh et ^ (l2013all u sed the photometry results 
from lCao et al.l (|2013ll . which are mag fainter 


than those we pres ent in this study based on 
lEldridge et al. ( 2015[l . While the predicted optical 
brightness is within the error bar of the observa¬ 
tion, the final mass of this model (^11 Mq) is 
too high to explain the typic al ejecta masses of 


SNe Ib/c (MB^pr-t.a. = 1-5 Mf^- Drout et'H] l201lL 
C and [ 201 ^ Taddia et al.l 201511 as well as the es¬ 


timated ejecta mass of iPTF13bvn (^1.9-2.3 Mq; 


‘^Systems where mass transfer from the primary star starts 
during the helium core contraction phase. 

^Systems where the first mass transfer starts during core 
hydrogen burning, followed by another mass transfer phase 
during helium core contraction. 


Bersten et al.ll20l4 iFremling et al.ll2014f) . 

Based on the SN ejecta mass and the opti¬ 
cal brightness of the progenitor candidate, a bi¬ 
nary progenitor with an initial mas s of 10-20 Mfn 


for iPTF13bvn has been sug gested l|Bersten et al 
2014 : Eldridee et al.l l2015li . Our stellar at¬ 
mospheric modeling of low-mass helium stars 
(Mflnai = 3-4 Mq) with an 0-type star compan¬ 
ion of ~ 20 - 30 M(n also supports th is scenario 
(Figures [5] and O. Bersten et al. ( 2014ll proposed 
the binary system composed of a 3.7 Mq helium 
star (Teff 16,000 K) and a 33.7 Mq hot com¬ 
panion (Teff 44,000 K) and predicted the SED 
of the proge nitor at the pre-SN s tate. Their model 
(Figure 5 of iBersten et al.l[2014 1 is fairly well con¬ 
sistent with the HST/AGS observations, but the 
assumption of a blackbody for the helium star 
may have uncertainty up to ^1 mag as discussed 
in Secion[5] Figure 0 compares a predicted spec¬ 
trum of one of our helium star binary models 
(M4w2-|-Ostar2) that reproduce best the obser¬ 
vations of the iPTF13bvn progenitor candidate, 
which are dereddened by the low extinction value 
{E{B — V) = 0.0715). In the composite spec¬ 
trum of the progenitor model, while the 0-type 
star dominantly contributes the flux in the optical 
wavelengths, emission lines expected from the he¬ 
lium star appear in the spectrum; this makes the 
helium star binary system to be distinguishable 
from 0-type stars with similar luminosity. 

7. Conclusions 

We have presented atmospheric models of rela¬ 
tively low-mass helium stars (3.0, 4.4, and 5.1 Mq) 
at the final evolutionary stage, which may repre¬ 
sent typical SN Ib/c progenitor s in binary systems . 


We confirm the prediction by Yoon et al. (20 I2I) 


that these binary progenitors can be significantly 
brighter than more massive WR progenitors in the 
optical, and the visual brightness becomes higher 
for a less massive helium star progenitor because 
of the more extended envelope. 

Their absolute magnitudes in the optical bands 
are comparable to those of 20 - 30 Mq 0-type stars 
{My = —5 ~ —6 mag). But unlike 0-type stars, 
their spectra could be marked by strong emis¬ 
sion lines if the mass-loss rate is sufficiently high. 
The envelopes of helium stars of 3 - 4 Mq un¬ 
dergo rapid expansion during the late evolutionary 
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Fig. 6.— Color-color diagram of (F435fF — 
F"555fF) vs. (F555VF—F814fF) of the helium star 
models without (filled diamonds) and with (filled 
diamonds with yellow borders) a 30 Mq 0-type 
star companion. Light and dark green filled stars 
are the colors of the iPTF13bvn progenitor can¬ 
didate dereddened by the low and high extinction 
values (see the caption of Figure 0, respectively. 



Fig. 7.— Predicted spectrum of the binary 
progenitor model composed of a 4.4 Mq he¬ 
lium star and a 30 Mq 0-type star companion 
(M4w2-|-Ostar2 model). Pink and orange lines 
present the model spectra of the helium star and 
0 -type star, respectively, and the black line is the 
composite spectrum of the two. The HST/ACS 
observ ations of the iPTF1 3bvn progenitor can- 
diate ( Eldridee et al. I l2015ll dereddene d by low- 


extin ction value of E{B — V) = 0.0715 (jCao et al 


201311 are compared with red squares. 


stages to byome a helium giant (77 ~ 10 — 5 0 Rq ; 
Yoon et al. I I 2 OIOII 2 OI 2 I: lEldridee et al.ll2015ll . and 


the resulting emission lines would be fairly narrow, 
compared to those of typical WR stars (Figures [2] 
and [3]). 

We compared our results with the observational 
properties of the progenitor candidate of the SN Ib 
iPTFlSbvn. We find that models with 3.0/4.4 Mq 
helium star plus a 20/30 Mq 0-type star com¬ 
panion give the best fit with the observation in 
term s of magnitude s and c olors, in good agreement 
with Bersten et al. ( 2014[l . but a 3.0 Mq helium 
star progenitor can also have optical magnitudes 
comparable to those of the iPTFlSbvn progenitor 
even without an 0-t ype star companion (see also 
Eldridee et al.l[2015i for a similar conclusion), de¬ 


pending on the degree of extinction to the source. 


Future observations might find the surviving O- 
type star companion of the iPTF13bvn progenitor. 
According to our prediction, its optical bright¬ 
ness should be lower by about AM = 0.2 — 0.9 
mag than that of the observed brightness of the 
iPTF13bvn progenitor, depending on the combi¬ 
nation of the helium star and 0-type star masses. 
This is because the contribution from the helium 


star progenitor must have disappeared. If no 
meaningful change in the optical magnitudes of 
the source is found in the future, it may indicate 
either that the iPTF13bvn progenitor was more 
massive than about 5 Mq while the source was the 
0-type companion star of 35 Mq, or that the 
source was not associated with iPTF13bvn at all. 
A caveat in this argument is that the luminosity of 
the companion star might have been significantly 
influenced by the interaction with the SN ejecta, 
and that there even exists the possibility that the 
source appears somewhat more luminous than in 
the pre-SN image, b ecause of the shock hea ting in 
the companion star ( Hirai fc Yamad£]l2015[) . 

However, as discussed above, a 3 Mq helium 
star alone without a luminous companion can also 
explain the optical brightness of the iPTFlSbvn 
progenitor in principle and non-detection of the 
surviving companion star would not necessarily 
imply a single star progenitor because the com¬ 
panion could_be_aH^nt_dwarf^ta£_or a compact 
object ( Eldridee et al. 2013 : Yoonll2015l) . With re¬ 
versal of the SN order in a Case A binary system, 
such a relatively low-mass helium star progenitor 


10 




























































produced in a binary system might not even ha ve 
any companion at the pre-SN stage (jPolslll994h . 

One important question still remains to be an¬ 
swered: if binary progenitors can be so bright 
in the optical as predicted by our models and if 
they represent typical SN Ib/c progenitors, why 
have they been so elusive in the previous searches? 
The detection limits were Mb — —4.4 mag and 
Mr ~ —4.92 mag in the searches of the progen- 
i tors of SN 2002ap a nd SN 2010br, respectively 
( Eldridee et al.ll201^ . They would have been de¬ 
tected, if their masses were about 3 Mq (see Fig¬ 
ure El), even without a luminous companion. How¬ 
ever, all of the other searches had more severe de¬ 
tection limits (Mbvri ^ —6.5 mag), and the pro¬ 
genitors would have been detected only with an 
0-type star companion of M > 20 — 30 Mq. The 
previous difficulty in detecting a SN Ib/c progen¬ 
itor might have resulted from the fact that most 
SN Ib/c progenitors in binary systems do not have 
a luminous 0-type star companion. It should also 
be noted that SN Ic progenitors would be gener¬ 
ally more difficult to detect than SN Ib progen¬ 
itors. This is because SN Ib progenitors would 
have a fairly massive helium envelope that results 
in large radii at the pre-SN stage (hence high opti¬ 
cal luminosities), while SN Ic progenitors that are 
helium-deficie nt would remain v ery hot and faint 
in the optical ( Yoon et al. 2012 1. SN 2002ap was 
a SN Ic, and therefore it is not surprising that its 
progenitor could not be detected even with such 
a deep search. This issue should be addressed 
more carefully with a binary population synthe¬ 
sis model that fully takes into account the evolu¬ 
tionary effects on the helium star structure at the 
final evolutionary stage and the recent observa- 
tional constraint s on the SN Ib/c ejecta masses (cf. 


tional constraint s on tne blN ib/c eiecta masses (ct 
Kochaned3l2009l: lYoon et ^ 2012t Eldridee et al 


20131). 
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Table 1 

Model Parameters 


Model 

M* 

(Mq) 

T* 

(^o) 

(K) 

Teff'’ 

(K) 

(^o) 

Rphot*’ 

(Rq) 

1 

0 

(km 8-1) 


V d 




V h 

M3wl 

3.0 

4.45e-|-04 

15057 

15060 

30.97 

30.97 

6.57e-07 

386.51 

1.0 

7.26e-02 

9.09e-01 

1.55e-04 

1.33e-02 

3.71e-04 

M3w2 

3.0 

4.45e-|-04 

15057 

14990 

30.97 

31.22 

l.OOe-05 

289.88 

1.5 

7.26e-02 

9.09e-01 

1.55e-04 

1.33e-02 

3.71e-04 

M4wl 

4.4 

9.05e-|-04 

28596 

28590 

12.26 

12.25 

1.50e-06 

452.57 

1.0 

l.lOe-01 

8.72e-01 

1.22e-04 

1.33e-02 

4.00e-04 

M4w2 

4.4 

9.05e4-04 

28596 

22630 

12.26 

19.55 

2.00e-05 

339.43 

1.5 

l.lOe-01 

8.72e-01 

1.22e-04 

1.33e-02 

4.00e-04 

M5wl 

5.1 

1.12e-|-05 

50587 

49520 

4.35 

4.54 

4.88e-06 

1332.76 

1.0 

7.35e-06 

9.81e-01 

2.46e-04 

1.32e-02 

2.94e-04 

M5w2 

5.1 

1.12e4-05 

50587 

39430 

4.35 

7.16 

2.00e-05 

999.57 

1.5 

7.35e-06 

9.81e-01 

2.46e-04 

1.32e-02 

2.94e-04 

Ostarl 

20.0 

7.05e-|-04 

29402 

29400 

10.23 

10.22 

1.50e-07 

1714.39 

1.0 

7.01e-01 

2.80e-01 

3.48e-03 

1.03e-03 

l.OOe-02 

Ostar2 

30.0 

1.99e-|-05 

32088 

32090 

15.00 

14.41 

6.00e-07 

1718.73 

1.0 

7.01e-01 

2.80e-01 

3.48e-03 

1.03e-03 

l.OOe-02 

Ostar3 

35.0 

3.72e-|-05 

28610 

28610 

24.84 

24.82 

1.13e-06 

1458.47 

1.0 

6.46e-01 

3.36e-01 

1.99e-03 

3.99e-03 

8.29e-03 


Note. —All the parameters except Tgff and i?« are taken from the stellar evolueionary models of [Yoon_et_^^ (201^ and used as the input parameters for the 
CMFGEN calculations. T’gff and are the outputs from the calculations. 

^From the stellar evolution code lYoon et alj f^UlCI) without correcting the optical depth effects from the wind. 

^Effective temperature and photospheric radius from the CMFGEN calculations defined at the Rosseland optical depth = 2/3. 

^From the standard /3-type velocity law. /3 = 1 and 1.5 for wl and w2 models, respectively, 
fraction of hydrogen, helium, carbon, nitrogen, and oxygen. 


Table 2 

Model Atoms Used in CMFGEN Calculation 


Species 

Super Levels 

Full Levels 

H I 

20 

30 

He I 

45 

69 

He II 

22 

30 

C II 

40 

92 

C III 

51 

84 

C IV 

59 

64 

N II 

45 

85 

N III 

41 

82 

N IV 

44 

76 

N V 

41 

49 

O II 

54 

123 

O III 

88 

170 

O IV 

38 

78 

O V 

32 

56 

O VI 

25 

31 

Si III 

33 

33 

Si IV 

22 

33 

Fe III 

104 

1433 

Fe IV 

74 

540 

Fe V 

50 

220 

Fe VI 

44 

433 

Fe VII 

29 

153 


Note. —The second column denotes the 
total number of super levels, which mean 
the energy levels that was used for the 
calculation. In the calculation, some lev¬ 
els having similar excitation energies are 
grouped into a single super level, to save 
computing time. In the case of H l, for ex¬ 
ample, 30 levels are considered in the cal¬ 
culation. For the first 15 levels, each su¬ 
per level corresponds to the actual energy 
level of H I, but the super levels from 16 
to 20 correspond to the grouped levels for 
16-18, 19-21, 22-24, 25-27, and 28-30, re¬ 
spectively. Therefore, the total number of 
super levels used for the computation is 20, 
while the total number of the considered en¬ 
ergy levels is 30, which is given in the third 
column. 
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Table 3 

Absolute magnitudes and bolometric corrections of the models and the blackbody 
APPROXIMATIONS IN THE HST/ACS F^SSW, F555W, AND F814-W FILTERS. 


Model 

Mf435W M.F555W 

Mf814W 

BCir'435W^ 

BCf555W^ 

BCf814W 


(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

M3wl 

-5.87 

-5.71 

-5.54 

-1.01 

-1.17 

-1.34 

M3w2 

-5.85 

-5.73 

-5.65 

-1.03 

-1.15 

-1.23 

M4wl 

-5.21 

-4.98 

-4.68 

-2.44 

-2.67 

-2.97 

M4w2 

-5.53 

-5.42 

-5.30 

-2.12 

-2.23 

-2.35 

M5wl 

-4.13 

-3.79 

-3.63 

-3.75 

-4.09 

-4.25 

M5w2 

-4.83 

-4.57 

-4.54 

-3.05 

-3.31 

-3.34 

Ostarl 

-4.86 

-4.61 

-4.28 

-2.52 

-2.77 

-3.10 

Ostar2 

-5.79 

-5.54 

-5.22 

-2.72 

-2.97 

-3.29 

Ostar3 

-6.75 

-6.52 

-6.22 

-2.44 

-2.67 

-2.97 

In binary systems (from composite spectra) 

M3wH-Ostarl 

-6.23 

-6.05 

-5.83 




M3w2+Ostarl 

-6.22 

-6.06 

-5.92 




M4wl+Ostarl 

-5.80 

-5.56 

-5.25 




M4w2+Ostarl 

-6.00 

-5.84 

-5.66 




M5wl+Ostarl 

-5.31 

-5.03 

-4.75 




M5w2+Ostarl 

-5.60 

-5.34 

-5.17 




M3wl+Ostar2 

-6.58 

-6.38 

-6.14 




M3w2+Ostar2 

-6.57 

-6.39 

-6.20 




M4wl+Ostar2 

-6.29 

-6.05 

-5.73 




M4w2+Ostar2 

-6.42 

-6.23 

-6.01 




M5wl+Ostar2 

-6.00 

-5.74 

-5.44 




M5w2+Ostar2 

-6.16 

-5.91 

-5.68 




M3wl+Ostar3 

-7.15 

-6.94 

-6.68 




M3w2+Ostar3 

-7.14 

-6.95 

-6.72 




M4wl+Ostar3 

-6.98 

-6.76 

-6.45 




M4w2+Ostar3 

-7.05 

-6.86 

-6.60 




M5wl+Ostar3 

-6.84 

-6.60 

-6.31 




M5w2+Ostar3 

-6.92 

-6.69 

-6.43 




Blackbody'^ 

BB_M3 

-5.91 

-5.85 

-5.81 

-0.97 

-1.03 

-1.07 

BB_M4 

-5.32 

-5.10 

-4.85 

-2.33 

-2.55 

-2.80 

BB_M5 

-4.01 

-3.72 

-3.38 

-3.87 

-4.16 

-4.50 


®'Based on the blackbody fluxes for the given temperature (T*) and luminosity of 3.0, 4.4, and 
5.1 Mq models in Table IT! 
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